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ABSTRACT 

We propose a new method for cosmological parameters extraction using the baryon 
acoustic oscillation scale as a standard ruler in deep galaxy surveys with photometric 
determination of redshifts. The method consists in a simple empirical parametric fit 
to the angular 2-point correlation function u){6). It is parametrized as a power law to 
describe the continuum plus a Gaussian to describe the BAO bump. The location of 
the Gaussian is used as the basis for the measurement of the sound horizon scale. This 
method, although simple, actually provides a robust estimation, since the inclusion of 
the power law and the use of the Gaussian removes the shifts which affect the local 
maximum. We discuss the effects of projection bias, non-linearities, redshift space 
distortions and photo-z precision, and apply our method to a mock catalog of the 
Dark Energy Survey, built upon a large N-body simulation provided by the MICE 
collaboration. We discuss the main systematic errors associated to our method and 
show that they are dominated by the photo-z uncertainty. 

Key words: data analysis - cosmological parameters - dark energy - large-scale 
structure of the universe 



1 INTRODUCTION 

Primordial perturbations generated acoustic waves in the 
photon-baryon fluid until decoupling (z ~ 1100). At this 
time the photons decouple from the baryons creating a high 
density region from the original source of perturbation, at a 
distance given by the sound horizon length. This high den- 
sity profile shows as a peak associated to the sound horizon 
scale in the galaxies spatial two-point statistics in the config- 
uration space, and as a series of oscillations in Fourier space, 
which can be used as a cosmological standard ruler. Re- 
cently, this feature has been applied to constrain cosmologi- 



cal parameters using several methods ( Eisenstein et al 
2006a|bl 



Hiitsi 



2007: 



Percival et al. 2007 



Padmanabhan et al 
Okumura et al. 2008 Gaztanaga et al. 2009; Sanchez 



2005 



et al.|2 009 ) . Those observations together with WMAP7 and 
Type la supernovae yield to the concordance cosmology: a 



spatially flat and late time accelerated universe (Komatsu 
et al.|2010 i. 



In order to achieve high accuracy on the cosmological 
parameters, several galaxy surveys will be done estimating 



redshifts from photometric data, like DES (The Dark En- 



ergy Survey Collaboration|2005[), P anSTARR S (|Kaiser et al 
2000 1), LSST ( jTyson et al.||2003| ) or PAU ( |Bem'tez et aT 
|2009 1. This higher accuracy will be possible thanks to a 
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larger volume and number of observed galaxies compared 
to previous spectroscopic surveys such as SDSS ( |York et al.| 
|2000[ > or 2dF ( |Colless et al]|2001[ ), even if these photomet- 
ric redshifts (photo-z) will have lower precision compared to 
their spectroscopic counterparts. The photo-z error depends 
mostly on the range of wavelengths covered by the filters 
in which the observations are done, the number of filters 
and their photometric error. As an example, the CFHTLS 
survey found, in the range 0.2 z ^ 1.5, a redshift error 
<t z ~ 0.03(1 + z). This error was obtained with observations 
in the optical and the near infrared ( jllbert et aLp)06| ). An 
error of this magnitude represents an uncertainty in the ra- 
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dial position of the galaxies which makes it impossible to 
infer the true 3-dimensional clustering pattern. Therefore, 
the analysis of angular statistics, like the two-point angular 
correlation function uj(9) and the angular power spectrum 
Ce, is required for these surveys. 

In this paper, a new method based on an empirical 
parametrization of cj(9) in redshift shells is proposed. The 
goal is to recover the angle corresponding to the BAO scale 
as a function of the redshift, and obtain the properties of 
the dark energy from its evolution. 

The method is designed to be used as a standard ruler 
only, i.e., we do not try to use the whole shape of the cor- 
relation function or the power spectrum. The reason is that 
although the strict standard ruler is less sensitive to the 
cosmological parameters, it seems also more robust against 
systematic uncertainties, since the only observation to be 
used is the position of the BAO peak in the two-point angu- 
lar correlation function, and nothing else. It is unclear which 
measurement will be more sensitive at the end of the day, 
including all the systematic errors, whether the full descrip- 
tion or the standard ruler ( Rassat et al. 2008 ; Sanchez et al 
20081. 



The main difficulty of an approach using the position of 
the BAO peak as a standard ruler is the extent at which the 
ruler remains "standard" . The position of the BAO peak is 
subject to several effects that shift its location and reduce its 
contrast. The most important among these is the projection 
offset due to the redshift bin width used to slice the data. 
This ultimately changes the to(9) shape for wide enough bin- 
ning leaving a shoulder shape with no local maximum. 

Our parametric approach is nonetheless able to measure 
the BAO scale even for wide redshift bins, and correct for 
the shift due to the projection effect, recovering the input 
cosmology. The proposed method has been tested against 
theoretical calculations and also applied to mock surveys 
from N-body simulations. Nevertheless, it should be useful 
to study BAO as a standard ruler in any photometric red- 
shift survey. 



2 ANGULAR CLUSTERING 

In complete analogy to the spatial correlation function £(r), 
the angular correlation function u>(9) is defined as the excess 
joint probability that two point sources (e.g. galaxies) are 
found in two solid angle elements d£l\ and dQ.2 with angular 
separation 9 compared to a homogeneous Poisson distribu- 
tion (Peebles 1980). From this definition, it is easy to find 



the relation between £(r) and uj{9), 



u>(9)= / dzi4>(zi) 
Jo 



dz 2 <j>(z 2 ) £(r; z) , 



(1) 



where cj>(z) is the redshift shell selection function normal- 
ized to unity and z = (zi + z%)l%. This relation assumes 
the time evolution of the real-space spatial correlation func- 
tion is small within the redshift shell under analysis, i.e., 
£(r; zi) ~ £(r; Z2) (for a relation without approximations 
Matsubara et al.|2004 i. The comoving distance between 



two point sources for a spatially flat cosmology is (for a gen- 
eralization to other spatial curvatures see again [Matsubara] 
et al.|2004| | 



V X{zi) 2 + X(z2) 2 - 2x(zi)x(z2) cos ( 



(2) 



The comoving radial distance x(z), assuming a constant 
dark energy equation of state parameter w Q is defined as 



dz 



Ho Jo V^(1 + ^) 3 + "a(1 + ^) 3 ( 1+ ^ 



(3) 



As usual, c is the speed of light, Ho is the present time 
Hubble parameter, f2 m is the matter density parameter and 
SI a = (1 — ^m) is the dark energy density parameter. With 
the relation between lj(8) and £(r) in hand, one needs a 
model for the two quantities that appear on Eq. the 
spatial correlation function, including the effects of non- 
linearities and the selection function <f>(z) that takes into 
account observational effects. 

The spatial correlation funcion is given by, 



3 dk 
2^2 



k 2 jo(kr)P(k;z) 



(4) 



where jo{x) is the spherical Bessel function of zeroth order 
and P(k) the power spectrum. 

In linear theory, valid at large distances, Pt{k) = 



b z D\z)P(k, 



0), with b a constant bias factor and 



D(z) = S(z)/S(z = 0), the growth factor of the density con- 
trast, S(z). Since we are only concerned with the BAO peak 
position and not the clustering amplitude, the bias b between 
luminous and dark matter will be set to unity throughout 
our analysis. We are assuming, then, that it does not evolve 
much in a redshift shell and that it is scale independent, see 
section [5.4.5| as well as neglecting peak shifts in the 3-d cor- 
relation of galaxies/halos with respect to that of dark matter 



in front of projection effects (Smith et al. 2008 Sanchez et al 



2008 1 . This is a very good approximation for the scales we 
are dealing with, as has been measured in [Crcsswcll fc Per-| 
cival (20091. A brief discussion of the systematic error that 



is induced by this assumption can be found in section [574] 

We introduce nonlinear matter clustering only through 
a damping of the BA O wiggles in the linear spectrum, i.e., 



P L e 



I Bharadwaj 1996 Crocce & Scocci- 



marro 2008; M atsubara[20 08 Eis enstein et al.|2007| |Crocce 
& Scoccimarro 20061. In doing this we discard the contri- 



bution to co(9) of the additive mode-coupling term to P(k), 



discussed in Crocce & Scoccimarro ( 2008 1 (see also Sanchez 
|et al.|2009| ), what is justified because it has a negligible im- 
pact in the angular position of the BAO peak in particular 
in comparison to projection effects ( Crocce et al.|2010 1 . The 
quantity a v characterizes the scale where non-linear effects 
become significant. It is related to the linear power spectrum 



<r v (z) = 



6tt 2 



dkP L {k;z) 



-1/2 



(5) 



In turn, the covariance matrix of oj(9), Covg 



1 Some of the cosmological models used in this analysis are non- 
flat and others have non-constant w. In these cases, the full cal- 
culation has been used. 
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(ui(8)uj(8')) , for a given survey can be estimated by 
(e.g. |Cabre et al.|2007| |Crocce et aL|2010[ ), 



Cov 9e / = 

^ 2(21 + l)P(cos(6>))P(cos(6>')) 



C(l) + 



N/AQ_ 



(6) 



where f s k y is the fraction of the sky covered by the survey 
and the ratio N/AQ is the number of galaxies per unit of 
solid angle in units of srad -1 . What we call statistical error 
along this analysis is computed using this expression and 
applied to u>(0). A detailed analysis of these expressions and 
their applicability is given in Crocce et al. (20101. 



3 THE SOUND HORIZON SCALE 

The strength of the standard ruler method lays in the po- 
tential to relate rather straightforwardly the acoustic peak 
position in the correlation function of galaxies to the sound 
horizon scale at decoupling. However, these two quanti- 
ties are not exactly equal (Sanchez et al.||2008| |Padman- 



abhan fc W hite 2009] |Seo et al.||2009| |2008| |Seo fc Eisen 



stein 



20071. Thus, in any standard ruler analysis it is cru- 



cial to distinguish between the following two angular scales, 
8b ao = r 's/x( z )> the angular scale corresponding to the 
sound horizon at the drag epoch and 8 P , the location of 
the BAO peak in the angular correlation function, defined 
as the local maximum around the expected angular scale. 
The peak position in the linear angular correlation func- 
tion 8 P only approaches the sound horizon scale Obao for 
infinitesimal shells, but a residual difference which ranges 
from 1 to 2 % is found ( |Sanchez et al.| |2008), depending on 
the cosmological parameters. 

Moreover, if the survey under analysis is performed with 
photo-z, there will be no such thing as an infinitesimal shell, 
due to the uncertainty in redshift associated to the photo-z. 
On the other hand, the impact of the redshift space distor- 
tions on oj(8) is large even for large angular scales, as has 
been pointed out by Nock et al. (20101, contributing also to 
these effects. 

A large photo-z error (e.g. o z ^ 0.03), will require a 
wide redshift shell. But it is well known that this fact in- 
duces two effects, which are shown in Figure [I] The first 
one is that the peak position shifts towards smaller angles, 
and the displacement is larger for larger widths. Second, the 
amplitude of the peak gets reduced until the local maximum 
disappears (Loverde et al. 2008; S impson et al.|[2009 ), and 
only a shoulder in uj(8) remains. As an order of magnitude, 
we have computed that for a z = 0.01(1 + z) there will be a 
well defined local maximum from mean redshift z > 0.6 and 
in the case <r z = 0.03(1 + z) it will be there only for z > 1.3. 

These results suggest we should measure 8b ao using an 
alternative method to the position of the local maximum of 
oj(8). Currently, the ratio of the power spectrum to a smooth 
fit is used, to avoid the effect of the tilting function, which 
changes the peak position ( Seo et aL"||2008 2009 Crocce & 



Scoccimarro 2008 Smith et al.|2008 l. Here we present an al- 
ternative method based on the angular correlation function. 
This will be the subject of the subsequent sections. 
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Figure 1. The impact of non-linearities and projection effects 
on the local maximum of the angular correlation function. The 
continuous line is the linear prediction and the dashed line the 
non-linear model from section [2] The figure assumes a redshift 
bin centered at z = 0.4 with the upper lines (blue) correspoding 
to an infinitesimal shell, the middle lines (green) to a shell width 
of Az = 0.05 and the lower ones (red) to Az = 0.1. Two effects 
are seen. First, the position of the peak moves towards lower 
angles when Az grows. Second, its amplitude reduces until the 
local maximum disappears, but a plateau survives instead. The 
dominant effect is the projection, while the impact of the non- 
linearities is much smaller. 



4 METHOD TO RECOVER 6» B ao 

Trying to overcome the problems induced by the photo-z, we 
propose a new method based on an empirical parametriza- 
tion of u>(6), which can be used even when there is no local 
maximum. The most important points of the method, the 
parametrization of u)(8) and the correction for the projection 
effect due to the width of the redshift bin are explained in 
detail in sections |4 . 1 1 and |4 . 2| respectively. The full recipe to 
obtain the BAO scale as a function of redshift is as follows: 

(i) Divide the full galaxy sample in redshift bins. 

(ii) Compute the angular two-point correlation function 
in each redshift bin. 

(iii) Parametrize the correlation function using the ex- 
pression: 



uj(8) =A + B8 1 + Ce 



-(8-B FIT ) 2 /2v 2 



and perform a fit to uj(8) with free parameters A, B, C, 7, 
8fit and a. 

(iv) The BAO scale is estimated using the parameter 
8 fit, and correcting it for the projection effect (see Eq. [8] 
below). The BAO scale as a function of the redshift is the 
only parameter needed to apply the standard ruler method. 
The cosmological interpretation of the other parameters is 
limited, since this is an empirical description, valid only in 
a neighborhood of the BAO peak, as is described in sec- 
tion gT] 

(v) Fit cosmological parameters to the evolution of the 
corrected 8bao with z. 
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h £l M fib f2fc wq W a n s 

0.70 025 0.044 O00 -1.00 00 0.95 
0.68 
0.72 

0.20 
0.30 

0.040 
0.048 

+0.01 
-0.01 

-0.90 
-1.10 

-0.1 
+0.1 

1.00 

Table 1. Summary of the 14 cosmological models used to test the 
method. Where empty, the fiducial values (first line) are assumed. 



4.1 Parametrization of lj(0) 

The correlation function u{&) is parametrized as in Eq. Q. 
A power law is used to describe the shape of the correlation 
function before and after the peak, and a Gaussian to de- 
scribe the BAO feature. The mean of the Gaussian is used to 
locate the BAO scale. The parameter A takes into account 
that the correlation function can go to negative values after 
the BAO peak. The parameters B and C describe the rel- 
ative weight of the two terms of the parametrization. The 
parameter 7 is the index of the power law, while 9 fit and 
a are the mean and the width of the Gaussian, respectively. 

We have tested the goodness of this parametrization in 
a redshift interval which ranges from 0.2 to 1.4, for a wide 
range of widths of the redshift bins (from to 0.2) and for 
14 cosmological models, which are listed in Table [E] The 
parametrization is good within 1%, i.e., the fit is good (the 
value of xVndof ranges from 0.98 to 1.01, and the prob- 
abilities of the fit from 0.6 to 0.9) when the error in each 
point of the correlation function is ~ 1% for all bin widths 
and cosmological models. This precision is much better than 



that expected in any realistic redshift survey, since it is 
much smaller than the cosmic variance. Several examples of 
the parametrization can be seen in Fig. [2] Other functional 
forms have been tested, but all of them are more complex 
and show no improvements on the description of the data. 

There is a region of stability to choose the starting and 
end points of the fit, which will be described in more detail in 
section [5.4.11 This region cannot be made arbitrarily large, 
since the description cannot be good for very low angles, 
due to non-linearities and changes in the power law index; 
neither for very large angles, since the correlation function 
changes its slope after the BAO peak. 

4.2 Correcting for Projection Effects 

The projection offset discussed in Sec. [3] is also present for 
the parametrization method. If we want to extract cosmol- 
ogy from this measurement, it is necessary to correct the 
fitted 6fit for the projection effect to recover 9bao- 

Obao = ol Ofit, (8) 

where a could in principle be a function of redshift, bin 
width and cosmology. However, if we want to do this correc- 
tion in an unbiased way, it must be cosmology independent. 
We have tested the method for the 14 different cosmologies 
summarized in Table[l] and for different redshift bin widths. 

For each of the 14 considered models, each of the con- 
sidered redshifts and each of the considered bin widths, we 
compute the angular correlation function. We use Eq.[l]and 
include non-linearities (see section [2j|. The calculation in- 
cludes also the effect of the galaxy distribution with the 
redshift, as descibed in section [53] The redshift ranges from 
0.2 to 1.4 and the bin width from to 0.2, as has been previ- 
ously stated. We then apply the fitting method to obtain the 
position of the BAO peak. Two main results are obtained: 

• The shift of Qfit with respect to Qbao has a universal 
shape, independent of the cosmology. Thus the correction 
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function a depends only on redshift and bin width, but not 
on the cosmological model, i.e., a = a(z, Az), within a pre- 
cision < 1%. 

• After applying this universal shift, the recovered BAO 
scale agrees with the true value for all cosmologies to a pre- 
cision of sC 0.75%. In particular, for infinitesimal bin widths, 
the method is able to compensate the displacement of the 
BAO bump due to the tilting function, correcting the 2% 



difference which has been observed in Sanchez et al. (2008). 



These two results are presented in Figure [3] where we 
plot the evolution of the shift, taking the angular scale cor- 
responding to the sound horizon scale as the reference, for 
5 different redshifts, for several bin widths and for the 14 
cosmological models. The spread of the results is constant 
with the redshift bin width, and comes from a possible resid- 
ual dependence on the cosmology together with the intrinsic 
limitations of the method (see section 5.3 on systematic er- 



rors). The average value of this shift is the correction applied 
to the fitted 8 fit- After this correction, the true value for 
the BAO scale is recovered for any bin width and cosmology. 
Note that the absolute value of the projection effect changes 
with the cosmological model, since both the BAO position 
and the radial distance change. However, the relative effect 
is cosmology independent, as shown in Figure [3] 



As has already been noted (Simpson et al. 20091, the 



projection effect is much more pronounced at low redshift. 
This must be taken into account to choose the optimum size 
of the redshift bins. The correction must be kept as small as 
possible, given the limitations imposed by the photo-z pre- 
cision, to introduce small systematic uncertainties. We are 
using a Gaussian form for the photo-z uncertainty, because 
if the requirements on photo-z measurement for a survey like 
DES arc fulfilled, the effects of a possible non-Gaussianity 
are kept small. If this is not the case, the evaluation of the 
true width of the bin must be refined. This is another reason 
to maintain the correction small. 

Note that the BAO scale is recovered both for linear 
and non-linear theory for infinitesinal bin width. This can 
be seen in Figure|4j where the residual of the fit with respect 
to the theoretical 9b ao as a function of redshift is shown for 
the fiducial cosmological model. The recovered values are 
well inside the 0.75% precision that we quote as systematic 
error, represented by the dotted lines. This happens for all 
the cosmological models, showing that the method is robust 
and able to recover the theoretical BAO scale. 



4.3 Redshift Space Distortions 

There are other effects which are important in order to have 
under control the parametric description of the angular cor- 
relation function. In particular, redshift space distortions 
must be taken into account in the analysis of the BAO scale, 
since they result in an anisotropic correlation function. We 
need to replace £(r) in Eq. [l]by f(<x, 7r) with iv — \'2 — Xi, 
the radial separation, and a = 2^1X2(1 — cos#), the per- 
pendicular separation. Galaxy pairs separated by large ra- 
dial distances infall into each other which means that they 
are measured as contributing to bins of smaller radial sepa- 
ration (in redshift space) than in real (or true) space. This 
results in larger correlations (more pairs) at intermediate 
scales (ir < 20 Mpc/h) and lower correlation at larger sepa- 
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Figure 3. Evolution of the shift of the measured 0fit with the 
width of the redshift bin for the 14 cosmological models consid- 
ered in this work. The sound horizon scale b AO f° r each cosmo- 
logical model is taken as the reference. Two results are obtained. 
First, the correct sound horizon scale is recovered for any cosmo- 
logical model within 0.75% for infinitesimal redshift shells. Sec- 
ond, the shape of the shift is universal, and the spread is constant 
with the redshift bin width. Note that the dispersion between 
models is much smaller (^ 1%) than the shift due to increasing 
binwidth. At low redshifts, z < 0.3, the shift can be ~ 15% for 
Ajz ^ 0.1, while at high redshifts, z > 0.5, the shift saturates at 
around 7%. 



rations. The redshift space correlation becomes negative for 
n > 40 Mpc/h while the real space correlation remains pos- 



itive until 7t > 130 Mpc/h (see Fig. Al in Gaztanaga et al. 
|2009[ ). Recall from Eq. [Ijhow the angular correlation is just 
a mean over the radial distance within the corresponding 
redshift bin. This means that for narrow redshift bins (of 
width comparable to it ~ 100 Mpc/h) the resulting angu- 
lar correlation can be quite different in real and redshift 



space (Fisher et al. 1995 Padmanabhan et al. 2007 Nock 
et al.pOlO l. As the total number of pairs is conserved, the 
integral over all separations is the same in real and redshift 
space, so that one expect to find the same angular correla- 
tion only for broad redshift bins, where the contribution of 
the boundary to the integral is negligible. 

The resulting predictions for redshift space agree very 
well with MICE simulations (more details can be found 
in Crocce et al.|p010 l. Fig. [5] shows some examples of this 
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Figure 4. Residual of the result of the parametric fit Ofit with 
respect to the the theoretical value of the BAO scale as a function 
of redshift for the fiducial model. The open dots show the results 
obtained from the fit to the linear correlation functions and the 
solid dots show the results on the non-linear correlation functions. 
In both cases the theoretical values are recovered, well inside the 
quoted systematic error of 0.75%, shown by the dotted lines. 



redshift space modeling for linear theory for two of the red- 
shift bins used in the analysis presented above. Symbols 
(with 2% nominal errors) correspond to the linear model 
predictions when we include the redshift space distortions 
as well as photo-z distortions. Note how both redshift space 
and photo-z distortions change quite dramatically the am- 
plitude of the resulting angular correlation. 

The lines in Fig. ISlshow the results of a fit of our 
parametrization (in Eq. [7| to these new predictions. The 
corresponding fit values of 8 fit are shown in the label. As 
expected, there is a small but systematic shift of Ofit to 
smaller angles in the case of photo-z distortions. This is be- 
cause the effetive redshift bin is wider due to the photo-z 
errors. This effect is taken into account by the correction 
shown in Fig. [3] (once we use the equivalent width of the 
true redshift distribution). Crocce et al. (2010 1 shows that 



the model we are using to predict the errors also works in 
redshift space. 

As can be seen in Fig. [5j the parametrization absorbs 
very well the effects of redshift space distortions and the 
value recovered for 6fit agree well in real and redshift 
space despite the large difference in amplitude. This is not 
too surprising as we have already shown above that this 
parametrization works for different cosmological models, dif- 
ferent redshifts and also different redshift widths. 



5 APPLICATION TO THE DARK ENERGY 
SURVEY 

5.1 The Dark Energy Survey 

The Dark Energy Survey (DES) F] is a next generation sky 
survey aimed directly at understanding the nature of the 



0.01 



3 



0.005 



0.015 



0.01 



3 



0.005 



z=0.49, Az=0.06 

• Redshift space true-z, d FtT = 4.495° 

• Real space true-z, Q RT = 4.495° 

• Redshift space photo-z, Q FIT = 4.366° 

• Real space photo-z, Q FIT = 4.369° 




9 (degrees) 



z=0.76, Az=0.076 

• Redshift space true-z, B Fj7 = 3.097° 

• Real space true-z, 9 nr = 3.099° 

• Redshift space photo-z, G F/r = 3.039° 

• Real space photo-z, 9 F/r = 3.041° 




http:/ /www. darkenergysurvey.org 
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Figure 5. Effect of z-distortions and photo-z. Points (with nom- 
inal 2% errors) show the linear theory prediction, while lines cor- 
respond to a fit to our parametrization in Eq. [7] Labels indicate 
the corresponding values of 9 fit m each case. 



dark energy. It will use 30% of the available time on the 
Blanco telescope at CTIO, during 5 years. The survey is 
designed to produce photometric redshifts in the range 0.2 < 
z < 1.4, and will cover 5000 square degrees in the southern 
hemisphere. The filters to be used are g, r, i, z, Y. 

DES has been designed to exploit mainly 4 methods 
in order to study the dark energy. These methods are the 
galaxy cluster counting and spatial distribution of galaxy 
clusters, the shifting of the galaxy spatial angular power 
spectra with redshift, weak gravitational lensing measure- 
ments on several redshift shells, and distances to ~ 2000 su- 
pernovae. Using these techniques, a 5%-15% precision mea- 
surement in the equation of state parameter w from each 
of our methods, and a 30% measurement in the variation 
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Figure 6. Photo-z precision used in this analysis (continuous 
line), as compare to used bins position and width (dots). 



chosen as a compromise between statistics, expected preci- 
sion in the photo-z determination and shift of the BAO peak 
due to projection effect. Then, the redshift bins used in this 
analysis follow the expected precision in the photo-z mea- 
surement, being wider for low and high redshifts and nar- 
rower for mid redshift. The analysis is performed in the full 
redshift range, from z = 0.2 to z = 1.4. With this prescrip- 
tion we have 14 bins, represented also as dots in Figure [6] 
The position of the dot in the horizontal axis represents the 
bin center, and the position in the vertical axis represents 
the bin width, which has been chosen to be close to the ex- 
pected 1 — cr error in photo-z measurement. Much wider bins 
make the sensitivity too small due to projection effect (the 
amplitude of the correlation function becomes very small 
and the shift of the peak becames very large), while much 
narrower bins make the correlations between redshift bins 
too large, complicating the analysis and also diminishing 
the sensitivity. 



of w with time are expected. Combined, they provide both 
stronger constraints and a check on systematic errors. 



5.2 N-body Simulation Features 

We have developed and tested our method to recover the 
angular scale of the sound horizon using a large N-body 
simulation capable of reproducing the geometry (e.g. area, 
density and depth) and specifications of DES r] 

The simulated data was kindly provided by the MICE 
project team, and consisted of a distribution of dark matter 
particles (galaxies, from now on) with the cosmological pa- 
rameters fixed to the fiducial model of Table Q] The radial 
profile matches the DES expectation 



dN/dz oc (z/0.5) 2 exp -(.z/0.5) 1 



(9) 



the simulation covers 1/8 of sky (~ 5000 square degrees) 
in the redshift range 0.2 < z < 1.4, containing 50 mil- 
lion galaxies. All these numbers roughly match the DES 
expectations for the main galaxies. These data was ob- 
tained from the comoving output at z = of one of the 
largest N-body simulations completed to date, with comov- 
ing size Ltox = 3072 /i" 1 Mpc and more than 8x 10 9 particles 
(m p — 2.3 x 10 11 Mo). More details about this run can 
be found inlFosalba et al.l (120081) andlCrocce et al.l (12009 1. 



The precision in the determination of the redshift for 
galaxies which has been used in this analysis is shown in 
Figure [6] This is the typical precision that is obtained using 
the current photo-z codes for DES (Ba nerji et al.|2 008). This 
photo-z has been introduced into the simulation, to produce 
a catalog of galaxies similar to that expected in DES, where 
the proposed procedure to recover Ob AO has been applied. 
The photo-z has been introduced by scattering each indi- 
vidual redshift according to the Gaussian representing the 
photo-z error of Figure [6] The redshift bin widths have been 



This catalog is publicly available at http://www.ice.cat/mice 



5.3 Results 

In order to numerically estimate the angular correlation 
function in a given redshift bin, we have used the Landy- 
Szalay estimator (Landy & Szalay 19931. The random cat- 
alogs contain as many galaxies as the simulation, since the 
shot noise error is negligible, given the size of the sample. 

The resulting correlation functions together with their 
parametrizations are presented in Fig. [jj The statistical er- 
ror in 9 fit is given as the error in the fit coming from the 
error in the correlation function, as explained in section [2] 
and taking into account the photo-z effect, as described in 
Eq. [12] below. 

The recovered values of Qbao after applying the correc- 
tion of Eq. [8] can be seen in Fig. [8] (top) as a function of the 
redshift. The correction is applied after evaluating the true 
width of the redshift bin, using Az tr ue = V^^ z photo, which 
corresponds to the width of a top-hat bin in true-z which 
gives the same amplitude in to(9) than the chosen photo-z 
bin. The points include also the systematic errors, described 
in section |5.4| The MICE cosmology is also shown as the 
solid line, and the best fit cosmology as the dashed line. 
To compare, the same analysis has been applied to the true 
redshifts catalog (i.e., without photo-z). The corresponding 
result is presented in the bottom panel of Figure |SJ together 
with the photo-z. The results obtained using the true red- 
shift are closer to the MICE cosmology, as expected. For true 
redshifts, the bin width is not corrected, and the systematic 
errors do not include any photo-z contribution. 

The true cosmology is recovered after applying the anal- 
ysis method, both for the photo-z results and for the true 
z results. The photo-z result is shown in Fig. [9] where the 
MICE cosmology is inside the l-cr contours. These contours 
have been obtained fitting the cosmology to the results pre- 
sented in Fig. [8] (top), with free parameters Q.m and w, and 
fixing all the other parameters to their MICE values. 

In order to perform this fit, we need to take into ac- 
count the correlation between bins induced by photo-z. To 
do this, we compute the migration matrix rtj, given by the 
probability that a given galaxy at bin i is measured at bin j 
due to the photo-z error. Once the migration matrix is com- 
puted (just counting the number of galaxies that migrate 
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from bin i to bin j in the simulation) , the calculation of the 
covariance matrix dj is straightforward. Since 



N 



3 = 1 



(10) 



where N, 



o 



is the number of observed galaxies in bin i, 

T 



the migration matrix and Nj is the true number of galaxies 
in bin j, then the observed and true correlation functions 
are related by ( Benjamin et al.|2010 |: 



u°{9) = J2 r l 



i.NP) 



(11) 



where wf '(9) is the observed correlation function in bin i 
at scale 9 and wf(9) is the true correlation function in bin 
i at scale 9, given by Eq. [j] Then, the covariance for the 
correlation functions can be written 



Cij = (Aw°(9)Aw°(8')} = (^40 



(iVj) 4 
(N°)i(N?) 



; Covg 



(12) 

where Covgg/ is given in Eq. [6] This expression neglects the 
intrinsic correlations between bins. Whether this is a good 
approximation depends mainly on the bin widths. In our 
case, adjacent bin centers are separated by more than 200 
Mpc//i, and we find that the intrinsic cross-correlations can 
be safely ignored. Thus, with photo-z errors, the observed 
correlation functions and covariances are given by Eqs. [Tl] 
and |12| respectively, rather than Eqs. [l] and [5] 

The resulting correlation coefficients, computed from 
Cij, for redshift bins are depicted in Figure [TO] The induced 
correlation extends to 3 bins. This correlation matrix is con- 
verted to the covariance of 9b ao angles, dj, including the 
corresponding errors. 

Therefore, we use the Cij to compute the \ 2 f° r the fit 
in the usual way: 



M,W) 



BAO ' 



7 BAO 



)C-j\9> E 



^BAo) 



(13) 



where 9 B ao ls the angular scale corresponding to the sound 
horizon for cosmological parameters f!jf , w at the redshift 
of bin i, and 9 B ao = a 9 fit is the measured angle in bin 
i. All the results include the systematic errors in the mea- 
surements of 9 bao, described more in detail in the next sec- 
tion. The statistical error is correlated using the computed 
covariance matrix. Some of the systematic errors (coming 
from photo-z and redshift space distortions) are also corre- 
lated, while the others are uncorrelated and are just added 
in cuadrature to the diagonal of the covariance matrix. All 
these numbers match the DES requirements. The total area 
of the survey, together with the magnitude limits, which fix 
the galaxy density, contribute to the errors on uj{9). More- 
over, the photo-z requirements fix the level of correlation 
among redshift bins. 

If we restrict to a one-dimensional analysis, and fix 
f!jn = 0.25, then we obtain for the equation of state of 
the dark energy w = —1.05 ± 0.14. The precision on this 
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Figure 8. Evolution of the measured Qbao with the redshift 
(top). Comparison of the results obtained using photo-z and the 
true-z (bottom). The true-z results are closer to the MICE cos- 
mology, as expected. 
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Figure 9. Allowed region for the cosmological parameters at 68 
(dashed-dotted line), 95 (dashed line) and 99 (solid line) % C. L. 
The true cosmology, given by the dot, is recovered within l-cr . 



result depends on the photo-z error, which is the dominant 
source of systematic uncertainties. If the photo-z error is de- 
creased, the precision in w correspondingly improves. For a 
more optimistic estimation of Az ~ 0.03(1 + z) it becomes 
Aw ~ 0.10. 
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Figure 10. Correlation matrix for redshift bins. The correlation 
extends to 3 bins, and is taken into account in the fit of the 
cosmological parameters. 



Systematic error 


^BAO Correlated between bins 


Parametrization 


1.0% 


No 


Photometric redshift 


5.0% 


Yes 


Redshift space distortions 


1.0% 


Yes 


Theory 


1.0% 


No 


Projection effect 


1.0% 


No 



Table 2. Estimation of various systematic errors. Some of them 
are correlated between redshift bins (photo-z and z-distorions er- 
rors). The others are not. 



5.4 Systematic Errors 

We have identified several sources of systematic errors for 
this methodology: 

• Uncertainties in the fitting procedure (method error). 

• Residual errors due to the correction in 8fit coming 
from a possible residual dependence on the cosmological 
model (cosmology independence). 

• Uncertainties in the theory coming from the implemen- 
tation of non-linearities (modeling error). 

• Effects of the redshift space distortions on the BAO 
scale (z-distortions error). 

• Uncertainty in the redshift determination (photo-z er- 
ror). 

A brief quantitative summary of these effects is given 
in Table [2j that follows from the different analysis discussed 
below. 



5.4.1 Method Error 

To compute the systematic error associated to the 
parametrization method, we have done some further analysis 
on theoretical angular correlation functions with the same 
bin widths and central redshifts as those used in the analysis 



z = 0.43 A z = 0.06 
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Figure 11. Bbao dispersion as a function of the starting and 
end point of the fitted region. Only two bins are shown, where 
the maximum of the effect is observed, at redshifts 0.43 (top) and 
1.015 (bottom). All the bins behave in a very similar way, but 
with smaller variations. The observed variation in the results is 
slightly smaller than 1%, and a conservative 1% systematic error 
is associated to the fit method. 

of the MICE simulation. The error associated to the method 
comes from the possible influence in the obtained 8 fit of 
the range of angles used to perform the fit. To evaluate the 
error, we have varied this range for all 14 z-bins we have, 
and performed the fit for each range. 

In the decision of the range to be fitted, we have to 
choose a starting point at angles smaller than the BAO peak 
and an end point after the peak. By varying this two angles 
we can study how much the result vary with this decision. 
Results can be seen in Fig. |11| where the obtained 6fit 
is shown for different starting points and end points of the 
fit, for low and high redshift. These results show the bins 
where the maximum differences are observed. In all cases, 
the uncertainty is always of the order of, but slightly below, 
1%. We have conservatively assigned a 1% systematic error 
associated to the method, since this error is subdominant. 



5.4-2 Cosmology Independence 

The dispersion of all the cosmological models around a com- 
mon value for the projection effect correction is small, ^ 1%, 
as can be seen in Figure [3] This dispersion propagates di- 
rectly to the corrected 8 fit as a source of uncertainty. Con- 
servatively, we have asigned a systematic error of 1% asso- 
ciated to the independence of the cosmology. 



5.4-3 Modelling Error 

Also the error due to the uncertainty in the theory (non- 
linearities at the scale of the BAO peak) has been computed 
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obtaining a global error of 1%, estimated in a conservative 
way as the difference between the Ofit measured using lin- 
ear and non-linear to(0), for the same redshift bins of the 
analysis. The difference for infinitesimal bins is presented in 
Fig-H 



5.4-4 Photo-z Error 

To compute the systematic error associated to photo-z mea- 
surement, we have analysed the simulation using the true 
redshifts in the same way we do for the photo-z. We can 
estimate the deviation in Ofit due to photo-z errors com- 
puting the difference between the result obtained using the 
catalog with true-z and the result obtained using the catalog 
with photo-z. In Fig. [15] (top) we can see the distribution of 
those differences. Photo-z uncertainties do not introduce a 
significant bias, but they do introduce a dispersion around 
the central value. This dispersion is around 5%, which we 
assign as the systematic error coming from the photomet- 
ric error. We are assuming here that this error is redshift 
independent, which given the scope of the current analysis 



is reasonable, as shown in Fig. 12 (bottom). The differences 



between the results of the analysis using photo-z and true z 
are shown as a function of the redshift, and there is no in- 
dication of any redshift dependence. This estimation could 
be further refined by using many simulations and repeating 
this study for each redshift bin as many times as simulations 
we use. This is out of the scope of this work, where we can 
only indicate the order of magnitude of the systematic er- 
rors, but should be done for a real survey, since this will be 
the dominant source of systematic uncertainty. 



5.4.5 Other Systematic Errors 

There are some other potential systematic errors, like those 
arising from the uncertainty in the description of the redshift 
space distortions, the gravitational lensing magnification or 
those mainly associated to the used galaxy sample, as the 
bias of the population. 

The uncertainties in the description of the redshift space 
distortions can induce a systematic error in the determi- 
nation of the sound horizon scale. To evaluate it, we have 
changed the parameters of the z-distortions within their un- 
certainties, finding a value of AObao ~ 0.5 — 1%, in good 
agrreement with the effect seen in Fig. [5] 

At the BAO scales it is a good approximation to con- 



sider that the bias is scale independent ( | Cress well fc Percival 
|2009||Crocce et al.|20io) , and given the width of the redshift 
bins, we can safely use a redshift independent bias inside 
each one. The redshift dependence of bias may appear in a 
different normalization of the correlation function between 
different redshift bins. Since we do not care about the overall 
normalization and only about the position of the peak, we 
are not sensitive to b(z). We have tested that the change in 
overal amplitude of the correlation function does not affect 
to our determination of the sound horizon scale, even if the 
error in the correlation function is changed. Moreover, we 
have tested the effect of a scale dependent bias, introducing 
artificially the effect in the correlation functions, using an 
approximated Q-model with the determination of parame- 
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Figure 12. Cumulated Difference between the measured Ofit m 
the photo-z catalog and the true z catalog for the 14 bins of the 
analysis (top). The dispersion is 5%, that we use as the systematic 
error. This error is taken as constant with redshift, since there is 
no clear dependence when the differences are plotted as a function 
of redshift (bottom). 



in the fit region of u(0) ranges from 1% to 6%, and in- 
duces a change in the determination of the sound horizon 
scale below 1%, well within statistical uncertainties. 

Finally, the gravitational lensing magnification can af- 
fect the uj(0) ( Loverde et al.| [2008) introducing correlations 
between bins at different redshifts. However, the effect of 
the magnification is concentrated at very low angles, and 
the magnitude of the induced correlation is much smaller 
than the correlation induced by photo-z, which has already 
been studied and used in this analysis. Therefore, the po- 
sition of the BAO peak is free of systematics coming from 
magnification. Consequently, all this effect is not considered 
here. 



ters of Cresswell & Percival (2009). The bias variation with 



6 CONCLUSIONS 

We have developed a new method to measure the BAO scale 
in the angular two-point correlation function. This method 
is adapted to photometric redshift surveys, where the infor- 
mation along the line of sight is lost due to the photo-z pre- 
cision, and only the angular information survives, although 
it can be applied to any survey. Two main results are found. 
First, the sound horizon scale can be recovered from the non- 
linear angular correlation functions to a precision ^0.75% 
applying the parametric fit described in the text, for any cos- 
mological model for infinitesimal redshift shells (i.e. for the 
3D correlation function). This is not totally surprising. It can 
be understood if we approximate the 3D linear correlation 
as a sum of a BAO gaussian and a baseline that is well de- 
scribed by a power-law. If we neglect mode coupling effects, 
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the non-linear correlation is given by a convolution with a 
gaussian damping, which leave unchanged the Gaussian po- 
sition. Under this assumption our parametric fit should be 
able to recover the true BAO positions without any bias. 
Second, the shift of the BAO peak due to projection effects 
has a universal shape, does not depend on the cosmological 
model, and only depends on the redshift and the redshift bin 
width. This can be used to correct the result obtained for 
wide photo-z bins and recover the true sound horizon scale. 
The method has been tested with a mock catalog built upon 
a large N-body simulation provided by the MICE collabora- 
tion, with characteristics similar to those expected in DES. 
The true cosmology is recovered within 1-a. The correlation 
between redshift bins and a preliminary evaluation of the 
systematic errors have been included in this study, and we 
find that the most important systematic error arises from 
the photo-z precision. The method is very promising and 
very robust against systematic uncertainties. 

Note that our analysis over the MICE simulations was 
in real, rather than in redshift space. Over a comoving out- 
put, rather than over an evolving distribution. Over dark 
matter particles, instead of galaxies. We believe that all this 
simplifications are not essential limitations to the method 
presented here, as we have shown that both the modeling 
and the error analysis are quite generic and work in real, 
redshift and photo-z space. If anything, redshift space dis- 
tortions (and possibly biasing) produce larger amplitudes, in 
which case our results are conservative in the sense that more 
realistic simulations (or observations) will produce smaller 
error bars than the ones presented here. 
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